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Abstract 

A  simple,  global  model  for  the  optimal  combination  of  ultracapacitor  and  battery  for  electrical  energy  storage  is  developed.  The  goal  for 
this  hybrid  storage  technology  is  to  reduce  the  system  life-cycle  cost  by  making  use  of  an  ultracapacitor’s  claimed  long  cycling-life  in  order  to 
supplement  relatively  cheap,  but  cycle-limited  battery  storage.  The  model  is  built  up  of  two  independent  sub-models  that  allow  flexibility  in  the 
relative  proportion  of  system  energy  storage.  An  analysis  performed  in  this  paper  indicates  that  ultracapacitor/battery  storage  systems  may  be  cost 
effective  for  high-cycle  applications. 
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1.  Introduction 

Renewable  energy  systems  that  make  use  of  fluctuating 
prime-movers  (e.g.  wind  power)  produce  power  that  also  fluctu¬ 
ates  on  several  time-scales.  For  a  wind  power  system,  time-scale 
fluctuations  range  from  the  slowest,  such  as  fluctuations  over  the 
course  of  several  years  (e.g.  El  Nino  effects),  to  seasonal  vari¬ 
ations  (e.g.  in  some  areas  the  wind  blows  stronger  and  more 
consistently  in  the  Winter),  to  short  time- scale  fluctuations  on 
the  order  of  seconds  up  to  an  hour.  Though  the  grid-connected 
load  also  varies;  in  many  cases,  it  is  desirable  to  smooth  out 
these  generation  fluctuations  to  some  degree. 

There  are  three  main  techniques  that  can  be  applied  in 
smoothing  out  wind  power  fluctuations  (1)  Aggregation  of  spa¬ 
tial  diversity  both  within  wind-farms  and  between  wind-farms 
(or  single  installations)  (2)  Aggregation  of  generation  diversity 
that  uses  other  means  of  generation  that  may  possess  differ¬ 
ent  fluctuation  characteristics  or  that  may  be  throttled  and  (3) 
Storage  of  surplus  energy  in  order  to  release  the  surplus  during 
times  of  power  generation  deficit.  Different  storage  technologies 
have  differing  characteristics  that  make  the  particular  technology 
more  or  less  suitable  for  reducing  the  variability  over  particu¬ 
lar  time- scales  and  for  particular  scales  of  energy  amounts.  For 
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instance,  pumped  hydro  storage  is  particularly  well- suited  for 
reducing  variability  in  the  hours-to-days  time-frame  for  very 
large  amounts  of  energy,  while  batteries  are  more  well  suited  for 
the  second  to  hour  time-frame  for  much  more  modest  amounts 
of  energy.  It  may  be  beneficial  to  combine  storage  technologies 
into  systems  that  allow  each  specific  technology  to  concentrate 
on  its  “best  fit”  of  time-scale  and  energy-scale,  thereby  creating 
a  storage  system  that  can  deal  more  effectively  with  variations 
from  power  sources.  Besides  power  source  smoothing,  a  well- 
matched  storage  system  could  also  be  useful  in  a  variety  of  other 
applications  (e.g.  power  smoothing  for  variable  loads). 

As  one  example  of  power  variability,  take  for  instance  one 
source  (namely  tower  shadow)  in  an  example  state-of-the-art 
electricity  producing  wind  turbine  with  a  nameplate  rating  of 
2  mW  and  an  overall  design  lifetime  of  approximately  20  years. 
The  turbine  rotor  will  rotate  at  approximately  15rpm  at  rated 
speed.  As  the  three  blades  pass  in  front  of  the  tower,  they  will  pass 
through  a  small  area  of  reduced  will  speed  directly  in  front  of  the 
tower.  This  tower  shadow  effect  reduces  the  wind  turbine  output 
by  up  to  approximately  1%  and  lasts  for  up  to  approximately  5° 
of  rotation.  Assuming  a  square  wave  for  the  power  pulse,  the  net 
effect  of  this  tower  shadowing  causes  a  0.75  Hz  56  ms  duration 
20  kW  power  dip  for  this  example  turbine.  As  the  mean  wind 
speed  fluctuates,  so  too  will  the  turbine  rpm,  produced  power, 
and  hence  the  duration,  frequency,  and  magnitude  of  the  power 
dip;  however,  if  the  turbine  is  placed  in  a  favorable  wind  regime 
it  will  be  operating  at  rated  speed  for  much  of  the  time.  As  a 
conservative  estimate,  let  the  wind  turbine  operate  at  rated  speed 
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Nomenclature 

al,a2,a3  coefficients  for  exponential  fit 
BBS  baseline  battery  system 

Cf  cycles  to  failure  for  battery 

DOD  depth  of  discharge 
RMSE  root  mean  squared  error 
2x(0.1U  +  0.9  BS)  combination  system  with  2  times 
energy  capacity  of  BBS,  utilizing  10%  Ucap 
Ucap  ultracapacitor 


for  10%  of  the  year  (i.e.  a  10%  cycling  duty  cycle),  this  translates 
into  roughly  2.4  million  cycles  per  year  of  this  example  power 
dip. 

Sources  of  power  variability  in  a  wind  turbine  are  not  limited 
only  to  tower  shadow;  wind  field  non-uniformity  across  the  rotor 
disk  (i.e.  turbulence  and  shear)  and  wind  field  non-uniformity 
over  time  (i.e.  wind  gusts),  as  well  as  switching  operations 
within  the  electrical  (and  sometimes  mechanical)  components 
and  actuator  lag  in  the  turbine  all  contribute  to  short-term  power 
fluctuations.  Given  the  wind  turbine  20  year  design  lifetime,  a 
very  large  number  of  cycles  (in  excess  of  106)  could  certainly  be 
achieved  by  fluctuations  with  cycling  frequencies  on  the  order 
of  0.016  Hz.  This  cycling  with  a  period  of  about  1  min  falls  in 
a  relatively  high-energy  portion  of  the  wind  spectrum  [1-3], 
and  so  a  high  cycling  duty  cycle  is  likely.  A  smoothing  system 
could  be  envisaged  that  would  reduce  the  magnitude  of  power 
fluctuations  resulting  in  a  more  desirable  constant  power  out¬ 
put  on  the  short-term  time- scale.  A  sophisticated  design  will 
take  all  sources  of  variability  into  account.  In  other  words,  a 
short-term  storage  system  should  be  designed  to  balance  the 
desires  to  smooth  power  over  a  range  of  fluctuations:  from  those 
induced  by  tower  shadow  (low  energy,  high  cycling  frequency) 
to  those  induced  by  wind  gusts  (moderate  energy,  low  cycling 
frequency),  perhaps  even  up  to  those  caused  by  slightly  longer- 
term  fluctuations  in  the  10-min  mean  wind  speed  (high  energy, 
very  low  cycling  frequency). 

This  paper  examines  the  possibility  of  some  optimal  combi¬ 
nation  of  battery  and  ultracapacitor  for  energy  storage  in  a  direct 
current  bus.  The  utility  of  the  lead  acid  battery  is  that  it  is  ubiq¬ 
uitously  available;  it  has  a  relatively  good  ratio  of  cost  to  useable 
energy,  and  it  is  familiar  for  use  in  renewable  energy  and  power 
systems  applications.  One  of  the  drawbacks  of  lead-acid  bat¬ 
teries  is  that  repeated  charging  and  discharging  leads  to  battery 
failure.  Due  to  (among  other  things)  repeated  cycling,  renewable 
energy  systems  can  be  particularly  troublesome  for  lead-acid 
battery  storage  [4].  Battery  construction  can  be  tailored  to  the 
application  so  that,  for  example,  with  thicker  lead  plates  battery 
longevity  may  be  improved  in  high  or  deep  cycle  applications. 
However,  it  is  unfortunate  that  though  batteries  are  rechargeable, 
using  them  in  this  manner  shortens  their  lifespan — in  some  cases 
dramatically. 

Like  batteries,  ultracapacitors  also  store  electrical  energy  via 
a  chemical  process;  however,  there  are  several  important  dif¬ 
ferences  between  the  two  storage  technologies.  Ultracapacitors 


Fig.  1.  Schematic  of  hypothetical  storage  integration. 

are  a  rather  new  technology  compared  with  lead- acid  batteries. 
Manufacturers  claim  that  cycling  (in  and  of  itself)  has  essen¬ 
tially  no  effect  on  the  lifetime  of  ultracapacitors.  Ultracapacitors 
can  discharge  at  much  higher  rates  than  batteries — though  only 
for  limited  durations.  Yet,  ultracapacitors  are  significantly  more 
expensive  than  batteries  on  a  cost  per  stored  energy  basis  and  do 
not  possess  the  high  volumetric  energy  density  of  lead-acid  bat¬ 
teries.  The  strengths  and  weaknesses  of  the  two  technologies 
should  fit  well  into  a  storage  system  where  the  ultracapac¬ 
itor  could  compensate  for  frequent,  short,  and  high  power 
disturbances,  while  the  battery  could  provide  compensation 
for  longer-term  less  frequent  events.  Fig.  1  shows  a  simple 
schematic  of  a  hypothetical  system  where  the  storage  has  been 
inserted  into  a  dc  bus  of  a  back-to-back  type  power  converter. 

2.  Methods 

The  overall  model  has  been  built  primarily  as  an  extension 
of  a  simplification  of  an  exponential  decay  and  Miner’s  damage 
rule  model  described  as  part  of  the  UMass  battery  model  [5]. 
This  double  exponential  decay  model  of  has  been  reduced  to  a 
single  exponential  decay  for  simplicity,  so  now  the  battery  life 
prediction  equation  for  cycles  to  failure  as  a  function  of  DOD 
uses  only  three  coefficients: 

CF  =  ai  +  a2ea3DOD  (1) 

Best-fit  coefficients  (for  a\,  02,  and  03)  for  the  BAE  OPzS  bat¬ 
tery  from  [5]  were  obtained  by  minimizing  the  RMSE  using 
the  Quasi-Newton  with  BFGS  Hessian  update  method.  Table  1 
shows  the  coefficient  values  obtained  and  95%  confidence  inter¬ 
vals  for  the  best-fit  coefficients.  Table  2  lists  the  datapoints  used 
in  creating  the  fit.  Fig.  2  shows  the  datapoints  (green  crosses 
-  DataPoints),  the  best  curve  fit  obtained  (blue  line  -  BF),  and 
also  the  curves  using  95%  confidence  interval  for  the  best- fit 
coefficients  (red  lines  -  95  Cl). 


Table  1 

Coefficient  Values  and  95%  Confidence  Intervals 


Coefficient 

Value 

95%  Cl 

a\ 

1.405  xlO3 

0.870  xlO3  to  1.941  xlO3 

a2 

1.341  xlO4 

0.947  xlO4  to  1.735  xlO4 

a3 

-7.312 

-9.957  to  -4.675 
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Table  2 

Data  Values  for  Example  Battery  [5] 


DOD 

CF 

0.1 

8000 

0.2 

4000 

0.3 

3200 

0.4 

2400 

0.5 

1900 

0.6 

1600 

0.7 

1350 

0.8 

1200 

This  example  battery  will  be  used  as  the  baseline  system 
in  the  following  analysis,  in  order  to  normalize  the  results  and 
allow  for  their  generalization.  Therefore,  DOD  will  be  fixed 
at  the  amount  for  the  baseline  system.  Fig.  3  shows  a  com¬ 
parison  of  the  effects  of  modifying  the  baseline  system  by,  for 
example,  doubling  the  size  (thereby  halving  the  effective  DOD) 
or  by  including  an  Ucap  as  a  percentage  of  the  system  total 
energy. 

In  the  figure,  the  baseline  battery  system,  BBS,  (blue  circled 
line)  possesses  a  shorter  cycling-life  than  a  battery  only  system 
that  is  twice  as  large,  2  x  BBS,  (red  circled  line),  a  system  that 
combines  a  Ucap  and  battery  in  a  ratio  of  0. 1-0.9  of  baseline 
energy,  1  x  (0.1U+  0.9BBS),  (blue  crossed  line),  or  a  combina¬ 
tion  system  of  0.2-1. 8  of  baseline  energy,  2  x  (0.1U  +  0.9BBs), 
(red  crossed  line).  Also  evident  in  the  figure  is  the  schema  for 
the  operation  of  the  combination  system:  discharges  of  less  than 
the  energy  level  of  the  Ucap  are  simply  assumed  to  not  affect 


Fig.  3.  Cycles  to  failure  of  systems. 


the  cycling-life  of  the  combination  system  e.g.  for  DODs  less 
than  0.1  for  the  1  x  (0.1U  +  0.9BBS)  system  or  for  DODs  less 
than  0.2  for  the  2  x  (0.1U  +  0.9BBS)  system. 

Costs  for  the  individual  components  of  the  model  are  sim¬ 
ply  taken  as  initial  capital  costs  and  are  linear  combinations 
of  current  prices  of  commercial  deep-cycle  lead-acid  batteries 
and  symmetric  ultracapacitors  and  are  approximately  5  x  10-5 
US$/Joule  and  1.5  xlO-2  US$/Joule,  respectively  [6,7].  The 
system  costs  are  also  normalized  with  respect  to  the  baseline 
system.  The  model  does  not  use  economic  methods  such  as 
net  present  value,  though  these  could  be  readily  added.  As  an 
example,  the  cost  for  the  2  x  BBS  system  is  2  times  that  of 
the  baseline,  while  the  1  x  (0.1U  +  0.9BBS)  system  would  cost 
approximately  35  times  the  baseline  system. 

There  are  two  main  design  variables  in  the  system  optimiza¬ 
tion:  ratio  of  Ucap  energy  to  system  energy  and  overall  system 
energy — relative  to  the  baseline  battery  only  system.  By  utiliz¬ 
ing  different  energy  ratios  and  system  sizes,  a  wide  range  of 
cycling  lifetimes  and  costs  can  be  explored  for  the  combination 
that  minimizes  the  lifetime  cost  of  the  system.  Both  the  battery 
and  the  Ucap  are  assumed  to  be  ideal  in  the  sense  that  only 
the  effects  of  DOD  on  cycling  life  are  being  studied  in  a  non¬ 
implementation  specific  design.  This  precludes  over  and  under 
voltage  operation,  partial  re-charging,  and  multi-cell  effects.  It 
also  does  not  take  into  account  any  efficiencies  or  reduced  oper¬ 
ational  abilities  for  either  battery  or  Ucap.  Finally,  while  the 
effect  of  cycling  on  the  battery  is  obviously  included,  the  Ucap 
is  assumed  to  be  unaffected  by  charge/discharge  cycles. 

3.  Results 

Before  discussing  the  results  of  the  system  modeling,  a  brief 
discussion  of  the  results  of  the  curve-fit  for  the  baseline  battery 
system  will  be  presented. 

The  curve-fit  problem  is  solved  by  minimizing  the  root  mean 
squared  error  between  the  candidate  functions  (i.e.  Eq.  (1)  with 
different  coefficient  values)  and  the  actual  datapoints.  Using  a 
Quasi-Newton  search  with  BFGS  update  method  in  the  Mat- 
lab  Optimization  Toolbox  [8]  required  36  main  loop  executions 
and  184  function  evaluations  from  a  start  of  [0  0  0]r.  Using  the 
steepest  descent  method  to  update  the  Hessian  did  not  result 
in  convergence  to  a  solution  within  an  acceptable  number  of 
iterations  (i.e.  greater  than  104).  Parallel  analysis  also  using  a 
least  squares  approach,  but  with  Levenberg-Marquardt  search 
utilizing  the  Matlab  Statistics  Toolbox  [8]  required  13  main  loop 
executions  from  a  start  of  [000]r;  resulted  in  the  same  coeffi¬ 
cient  estimates;  and  yielded,  using  the  linearized  estimate  of  the 
coefficient  covariance,  95%  confidence  intervals  (see  Table  1) 
for  the  predicted  coefficient  values.  Several  starts  from  different 
initial  conditions  resulted  in  the  same  optimal  final  combina¬ 
tion  of  coefficients,  presumably  this  is  predicted  to  some  degree 
by  the  fact  that  there  are  eight  datapoints  (equally  spaced  over 
the  v-axis)  and  only  three  coefficients  resulting  in  an  overcon¬ 
strained  problem  that  should  have  only  one  solution.  In  other 
words,  given  only  three  datapoints,  it  would  be  possible  to  fit 
a  unique  curve  that  would  possess  an  RMSE  of  zero,  since  this 
equation  has  an  additional  five  datapoints,  the  problem  becomes 
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Fig.  4.  Curve  fit  coefficients  and  log(RMSE). 


minimizing  the  RMSE.  The  objective  function  to  be  minimized 
is: 


f{a\,  a2,  <23) 


n 

-  (ai  +  a2efl3DOD')]2  (2) 

i=  1 


where  i  ranges  from  1  to  8  (the  number  of  datapoints).  At  an 
extremum  (e.g.  the  minimum  RMSE)  the  gradient  of  /  with 
respect  to  a\a2,  and  <23  should  vanish.  In  order  for  the  extremum 
to  be  confirmed  as  a  minimum,  the  Hessian  of/ at  (at  least)  the 
extremum  should  be  positive  definite.  Of  course,  there  is  no 
guarantee  that  this  minimum  is  anything  but  a  local  minimum 
unless  the  Hessian  can  be  demonstrated  to  be  positive  definite 
over  the  entire  space  of  a\,  <22,  and  <23.  The  definition  of  RMSE 
guarantees  a  positive  definite  Hessian  over  the  entire  space  for 
real  values — due  to  the  squared  term.  This  means  that  the  solu¬ 
tion  found  is  a  global  minimum  for  the  RMSE.  Fig.  4  shows 
a  three-dimensional  volume  plot  where  the  color  of  the  region 
indicates  the  natural  logarithm  of  the  RMSE  value  (blue  indi¬ 
cates  smaller  RMSEs  while  red  indicates  larger  RMSEs).  The 
log  function  has  been  used  in  the  plot  in  order  to  enhance  the 
contrast.  Also  shown  in  the  figure  is  the  optimization  trajectory 
(the  black  dotted  line).  As  expected,  the  trajectory  points  toward 
the  smallest  RMSE  value. 

Because  the  coefficients  are  constants,  the  sensitivity  of  the 
objective  function  (with  respect  to  the  individual  datapoints)  is 
only  a  function  of  the  partial  derivative  with  respect  to  each  data- 
point.  When  one  differentiates  the  objective  function,  in  essence, 
one  recovers  the  original  equation  for  battery  life,  i.e.: 


d 

dDOD; 


/ 


=  a2ajea3DODiqqi 


d 

dCF,i 


/  =  —qqiC-pj 


(3) 

(4) 


where 


qqi  =  -Cf,/  +  a\  +  tf2e^3DOD/ 


(5) 


Assuming  that  the  terms  common  to  all  variables  are  factored 
out. 

Because  the  sign  of  the  exponent  of  the  exponential  term 
(i.e.  <22)  is  negative,  the  sensitivity  analysis  of/ with  respect  to 
DOD  shows  that  the  objective  function  seems  to  more  heavily 
weight  datapoints  at  smaller  values  along  the  abscissa  (i.e.  data 


Fig.  5.  Relative  cost  of  combination  system  compared  to  life  comparable  battery 
only  system:  battery  only  system  (blue  Line),  10%  Ucap  system  (green  circles), 
20%  Ucap  system  (red  crosses),  50%  Ucap  system  (purple  triangles). 

lying  closer  to  the  origin).  Fortunately,  the  fitted  curve  in  this 
case  not  only  possesses  a  minimal  RMSE,  but  it  also  seems  to 
qualitatively  fit  the  data  very  well. 

Fig.  5  shows  a  simple  comparison  between  the  costs  of  a  Ucap 
plus  battery  system  and  a  battery  only  system  that  is  sized  to 
provide  comparable  cycle  life.  In  order  to  improve  legibility,  the 
figure  is  normalized  such  that  only  DOD’s  above  the  Ucap  ratio 
for  each  system  are  considered  (i.e.  the  curve  for  the  10%  Ucap 
system  has  been  shifted  by  —0.1,  and  so  on).  As  can  be  seen  in 
the  figure,  as  the  DOD  rises  above  the  capacitor  size,  the  relative 
costs  of  the  system  increase  rapidly.  However,  there  is  a  region 
(circled  in  red)  where  the  capacitor/battery  combination  should 
cost  less  for  a  similar  number  of  cycles  at  particular  depths  of 
discharge. 

A  more  complicated  picture  emerges  when  the  two  design 
variables  (Ucap/System  energy  ratio  and  System  size)  are  indeed 
allowed  to  vary  over  the  solution  space.  Figs.  6  and  7  illustrate 
the  system  cost  over  the  design  space,  and  are  three-dimensional 
views  of  function  values  for  DODs  of  0.9  and  0.1  respectively; 
both  the  height  along  the  z-axis  and  color  spectrum  indicates 
the  function  value  (i.e.  tall  red  indicates  systems  that  cost  more 
while  low  blue  indicates  systems  that  cost  less). 

Also  visible  in  the  figures  are  optimization  trajectory 
endpoints  (i.e.  optimal  points).  As  might  be  suspected,  the 
non-convex  nature  of  the  solution  space  proves  a  challenge 
for  gradient-based  optimization.  In  the  figures,  one  trajectory 
endpoint  (green  box)  possesses  a  initial  values  that  cause  the 
search  to  proceed  to  a  local  minimum,  while  the  other  trajectory 


Fig.  6.  System  cost  for  0.9  DOD. 
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Fig.  7.  System  cost  for  0.1  DOD. 
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Fig.  8.  Relative  cost  to  Ucap  failure,  0.5  DOD,  BF  coefficients. 

endpoint  (red  circle)  possesses  an  initial  condition  that  causes 
the  solution  to  proceed  to  the  global  optimum.  The  difference 
between  the  global  and  local  optimum  is  that  due  to  the  nature 
of  the  model,  the  cost  of  a  system  that  possesses  a  Ucap  that  is 
sized  in  order  to  match  the  DOD  approaches  zero  as  the  num¬ 
ber  of  cycles  (i.e.  cycle-lifetime)  tends  towards  larger  and  larger 
values;  while  the  baseline  system  cost  remains  at  unity.  This  is 
somewhat  an  artifact  of  the  model,  so  another  analysis  has  been 
performed  with  an  upper  limit  on  the  cycle-lifetime  of  the  Ucap. 

In  this  revised  analysis,  a  lifetime  limit  of  one  million  cycles 
has  been  imposed  on  the  Ucap  component,  and  batteries  in  the 
systems  are  replaced  once  their  lifetimes  have  expired.  In  effect 
the  question  asked  in  this  case  is  “How  much  does  it  cost  for 
this  system  to  reach  one  million  cycles?”  Figs.  8-10  present 
the  three-dimensional  results  from  the  case  where  the  Ucap  size 
is  allowed  to  vary  while  the  DOD  (again  using  the  baseline 
DOD)  is  fixed  at  0.5.  In  the  figures  the  black  region  is  the  region 


Fig.  10.  Relative  cost  to  Ucap  failure,  0.5  DOD,  HC  coefficients. 


where  the  cost  of  a  combination  system  would  be  equal  to  a 
comparable  battery  only  system;  therefore  systems  that  have 
costs  that  are  below  the  black  plane  are  more  cost-effective  than 
battery  only  systems  (with  exhausted  battery  replacement  in  both 
cases). 

Fig.  8  shows  results  using  the  best-fit  coefficient  values  - 
hence  the  caption’s  . . .  BF  Coefficients’ .  Figs.  9  and  10  use  coef¬ 
ficient  values  from  the  confidence  interval  limits.  Fig.  9  shows 
that  the  region  of  cost-effectiveness  for  a  Ucap/battery  system 
has  grown  (as  expected)  due  to  the  use  of  the  confidence  inter¬ 
val  coefficients  that  result  in  a  reduced  cycling  capability  for 
the  battery  portion  of  the  system  -  the  LC  (low  cycle)  coeffi¬ 
cients.  Fig.  10  also  shows  expected  results  in  that  the  region  of 
cost-effectiveness  has  shrunk  due  to  the  use  of  the  confidence 
interval  coefficients  that  result  in  increased  cycling  capability 
for  the  battery  portion  of  the  system  -  the  HC  (high  cycle)  coef¬ 
ficients.  Even  using  the  HC  coefficients,  there  is  still  a  region 
where  the  Ucap/battery  system  should  be  more  cost-effective 
than  a  strictly  battery -based  system. 

Neither  the  minimums  (with  Ucap  to  system  ratio  =0.5  and 
system  size  =1)  nor  the  fact  that  a  gradient-based  optimization 
algorithm  works  well  (for  the  convex  solution  space)  is  surpris¬ 
ing.  The  figures  are  representative  of  the  findings  for  all  other 
DODs. 

What  is  more  surprising,  however,  are  the  results  presented 
in  Figs.  11-14  .  In  this  analysis,  the  Ucap  to  system  energy 
is  fixed  at  given  values  (in  the  case  of  Figs.  11  and  12 
Ucap/System  energy  ratio  =  0.1)  and  the  DOD  is  allowed  to 
vary.  Figs.  13  and  14  show  the  results  of  a  similar  analysis  for 
energy  ratio  =  0.5.  The  analysis  has  been  performed  using  the 


Fig.  9.  Relative  cost  to  Ucap  failure,  0.5  DOD,  LC  coefficients.  Fig.  11.  Relative  cost  to  Ucap  failure,  0.1  Ucap  energy  ratio,  BF  coefficients. 
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Fig.  12.  Relative  cost  to  Ucap  failure,  0.1  Ucap  energy  ratio,  HC  coefficients. 
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Fig.  13.  Relative  cost  to  Ucap  failure,  0.5  Ucap  energy  ratio,  BF  coefficients. 


Fig.  15.  Relative  cost  to  Ucap  failure,  0.5  Ucap  energy  ratio,  projected  Ucap 
cost,  BF  coefficients. 
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Fig.  14.  Relative  cost  to  Ucap  failure,  0.5  Ucap  energy  ratio,  HC  coefficients. 


best- fit  coefficients  (Figs  11  and  13)  and  the  more  limiting  (in 
terms  of  cost-effectiveness)  HC  coefficients  (Figs  12  and  14). 

As  can  be  seen  in  each  of  the  figures,  there  is  a  region  (in 
black)  where  the  cost  of  a  relatively  small  amount  of  Ucap  stor¬ 
age  can  offset  the  cost  of  the  batteries  required  (and  therefore 
provide  a  lower  overall  system  cost)  over  a  wide  range  of  DODs. 
This  region  changes  shape  (though  seems  to  remain  convex)  as 
the  Ucap  to  system  energy  ratio  changes. 


4.  Discussion 


Ultracapacitor  cost  is  projected  to  decrease  [7]  and  the  effect 
of  cost  reductions  will  be  to  increase  the  size  of  the  region  of 
the  below  baseline  system  cost.  Fig.  15  shows  the  results  of 
the  high-cycle  system  cost  analysis  using  the  projected  cost  of 
ultracapacitors  in  a  system  identical  to  that  used  for  Fig.  13; 


Fig.  16.  Relative  cost  to  Ucap  failure,  0.5  Ucap  energy  ratio,  projected  Ucap 
cost,  HC  coefficients. 

the  systems  used  in  the  analysis  for  Figs.  14  and  16  are  also 
identical. 

When  Figs.  15  and  16  is  compared  with  Figs.  13  and  14,  it  can 
be  seen  that  the  area  of  applicability  has  significantly  increased. 
It  is  not  immediately  visible  in  the  figures  that  the  overall  system 
cost  compared  to  the  results  shown  in  Figs.  13  and  14  (i.e.  the 
values  on  the  z-axis)  has  also  decreased  substantially. 

5.  Conclusion 

The  results  of  the  fixed  Ucap-to-system  energy-ratio  with 
varying  depths  of  discharge  are  very  encouraging,  in  that  they 
imply  that  a  small  amount  of  Ucap  storage  may  reduce  over¬ 
all  system  lifetime  cost  (to  a  large  number  of  cycles)  over  a 
wide  range  of  depths  of  discharge.  This  ability  to  benefit  a 
range  of  DODs  will  be  important  in  applications  where  the  DOD 
can  vary  significantly,  for  example  short-term  renewable  power 
smoothing. 

High  cycling  applications  are  present  in  many  areas  of  the 
electrical  power  system.  The  utility  of  this  simple  model  is  to 
show  that  a  small  amount  of  Ucap  storage  may  well  reduce 
battery-based  smoothing  system  overall  lifetime  cost.  This  pro¬ 
vides  the  motivation  to  investigate  further  including  such  issues 
as  more  sophisticated  modeling  (for  each  component,  the  system 
as  a  whole,  and  the  loading  conditions  to  which  such  a  system 
will  be  subject)  and  methods  of  implementation  of  such  a  hybrid 
storage  system. 
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